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We report on an extensive study of the influence of spin anisotropy on spin glass aging dynamics. 
New temperature cycle experiments allow us to compare quantitatively the memory effect in four 
Heisenberg spin glasses with various degrees of random anisotropy and one Ising spin glass. The 
sharpness of the memory effect appears to decrease continuously with the spin anisotropy. Besides, 
the spin glass coherence length is determined by magnetic field change experiments for the first time 
in the Ising sample. For three representative samples, from Heisenberg to Ising spin glasses, we 
can consistently account for both sets of experiments (temperature cycle and magnetic field change) 
using a single expression for the growth of the coherence length with time. 

PACS numbers: 75.10.Nr,75.50.Lk,75.40.Gb 



Real spin glasses (SGs), which are both disordered and 
frustrated magnetic systems, are always found to be out 
of equilibrium on the laboratory time scale. When cooled 
below its freezing temperature T g , a SG starts to age; 
from a completely random initial configuration, spin-spin 
correlations develop and progressively 'rigidify' the spin 
network. If the SG is further cooled, the aging process 
restarts (rejuvenation effect) but the previous aging is 
not lost. When reheated, the SG retrieves its former 
state (memory effect) and resumes its former evolution. 
These memory and rejuvenation effects suggest a descrip- 
tion in terms of a hierarchical organization of metastable 
states as a function of the temperature 0, although an 
alternative description in terms of temperature chaos has 
also been proposed 0. Recently, a comparison between 
an Heisenberg and an Ising SG showed that the slow 
dynamic properties are quantitatively different in these 
systems |3|, the memory effect being more pronounced 
in the Heisenberg case. How are the spin anisotropy and 
the aging properties related? Is there a systematic effect 
of the anisotropy? Answering these questions would shed 
a new light on the longstanding and topical problem of 
the nature of the Heisenberg SG phase |4|. The relevance 
of anisotropy for this problem is supported by the recent 
finding that the critical exponents vary continuously with 
the SG anisotropy 0- 

We report in this letter on two classes of experiments. 
First, we have performed thermo-remanent magnetiza- 
tion (TRM) measurements with negative temperature cy- 
cles during aging in one Ising and several Heisenberg SGs 
in order to quantify the influence, at temperature T, of 
aging at a slightly smaller temperature T — AT. There 
is indeed a small cumulative effect from one temperature 
to the other which limits the sharpness in temperature 
of the memory effect. It is this cumulative effect that 
we measure here as a function of the spin anisotropy. 
Secondly, we have performed new magnetic field change 
experiments in the Ising sample. This protocol enables 
us to estimate the coherence length that grows during 
the aging process. Both sets of experiments are success- 



fully analyzed within the framework proposed in Ref. |6j 
(building upon ideas developed in 0, see also Q), which 
proposes a real space picture for the hierarchical organi- 
zation of the SG states. In this model, at a given tem- 
perature and time scale, aging processes are associated 
to a specific length scale. Larger length scales are frozen 
while smaller length scales are fully equilibrated. Upon 
reducing the temperature and as a result of the change 
in the Boltzmann weights, the small length scales which 
were equilibrated at higher temperature have to evolve 
towards a new equilibrium state (rejuvenation). On the 
other hand, the structure at large length scales can no 
longer evolve (memory). In this picture, the sharpness of 
the memory effect is directly related to the rapid varia- 
tion of the active length scale with temperature. A cen- 
tral concern of the present paper is to establish exper- 
imentally the time and temperature dependence of this 
active length scale (see Eq. 0] below.) 

The temperature cycle experiments were performed in 
a monocrystalline Feo.sMno.sTiOa Ising |9j sample #1 
(T g = 20.7K) and Heisenberg SGs with decreasing ran- 
dom anisotropy arising from Dzyaloshinsky-Moriya in- 
teractions: an amorphous alloy (Feo.iNio.g^sPieBgALj 
#2 (T g = 13AK), a diluted magnetic alloy Au:Fe 8% #3 
(T g = 23.9K), an insulating thiospinel CdCri. 7L10.3S4 
#4 (T g = 16.7K) and we reproduce the data on the 
canonical Ag:Mn 2 . 7% SG #5 (T g = 10.4K") from Ref. [ljj. 
The relative anisotropy constants (K /T g ) / '(K /T g ) Ag-.Mm 
measured by torque experiments |5|], are respectively 
16.5, 8.25, 5 and 1 for samples #2, #3, #4 and #5. 
An isothermal TRM procedure allows one to define a set 
of reference curves that depend on the waiting time, and 
to which we will compare the relaxation curves obtained 
after more complicated histories. In our series of temper- 
ature cycle experiments, the sample is first cooled to a 
temperature T < T g under a small magnetic field (a few 
Oe). After waiting a time t\ = 500s, the sample is further 
cooled to a slightly smaller temperature T — AT during 
<2 = 9000s and then re-heated to T for another short 
time £3 = ti (see inset in Fig. [3 ) ■ The additional waiting 
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FIG. 1: Effective waiting times deduced from the tempera- 
ture cycle experiments at T ~ 0.72T 9 (a, filled symbols) and 
T ~ 0.85T 9 (b, open symbols). The SGs are labelled from #1 
to #5 by order of decreasing anisotropy. 



times t\ and t% permit to control any possible thermal- 
ization delays. The magnetic field is then switched off 
and the TRM is recorded. For each 'cycled' TRM curve, 
an effective waiting time t 2 ^ is defined such as to super- 
impose the obtained TRM with one of the purely isother- 
mal TRM recorded at temperature T, with waiting time 
t w = ti + t L 2 ff (T, AT) + t 3 . For AT = 0, the effective 
waiting time equals the actual waiting time but becomes 
smaller for AT > 0: t%, ff (T,AT) < t 2 characterizes the 
effect at T of aging at T- AT. The results for the differ- 
ent SGs and for two different temperatures, T/T g ~ 0.72 
and T/Tg ~ 0.85, are plotted in Fig. ^ (for sample #5, a 
similar procedure, with ti = £3 = 0, was applied [f0()- 

The slope of t e 2 ^ /t 2 vs AT/T reflects the cumulative 
effect of aging from one temperature to the other. Hence, 
it quantifies the memory effect; the steeper the curve, the 
sharper the memory effect. For each of the temperatures 
T/Tg, the influence at T of aging at T — AT is found 
to increase with increasing spin anisotropy. Note that, 
at T/T g — 0.72, the measured for sample #3 are 
slightly smaller than those of the less anisotropic sample 
#4. However, these two samples are very similar as far 
as their anisotropy is concerned since their anisotropy 
constants differ only by 50% whereas from sample #1 to 
#5 they spread over more than one order of magnitude. 

As can be seen in Fig. [3] below, for a given AT/T, 
the effective time decreases with increasing temperature 
for all the Heisenberg-like samples (samples #4 and #5). 
By contrast, in the fsing sample #1, the influence at T 
of a T — AT aging is slightly stronger at the higher tem- 
perature in the Ising case. Previous ac-measurements 
on the fsing sample 3] demonstrated this singular be- 
haviour even more clearly. At the higher temperature 
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FIG. 2: Schematic representation in terms of SG correla- 
tion length of the negative temperature cycling experiment 
depicted in the inset. 



T/Tg = 0.84, the dynamics was in fact found to be al- 
most critical (nearly no AT dependence). 

To interpret these experiments in a real space pic- 
ture @, we introduce a typical length scale i(t w ,T) over 
which spin-spin correlations establish at temperature T 
and waiting time t w . From general arguments, this co- 
herence length I should be an increasing function of T 
and t w . For small AT, the experimental TRM curves 
obtained after a temperature cycle superimpose onto ref- 
erence isothermal curves. This implies that the SG has 
reached the same state after a t\ + t% isothermal ag- 
ing and the t\ + t 2 experimental temperature cycle. We 
can then sketch the temperature cycle experiment as in 
Fig. [3 where t e 2 ^ is defined by: 



£(tl+t 2 ,T-AT) 



£(t 1+ t e 2 ff ,T) 



(1) 



where t\ is given by l{h,T) = £{t\,T - AT). Note that 
this picture which retains only the cumulative effect of 
aging with temperature (no rejuvenation effect) is only 
valid for small temperature cycle AT. The slope of the 
experimental curve t e 2 fi ' /t 2 vs (AT/T), which quantifies 
the sharpness of the memory effect, reads from Eq. Q 
and for AT/T < I and h < t 2 : 



f eff 



/ dl/dlnT 
\d£/d\nt w 



AT 
~T~' 



(2) 



As announced, a sharp memory effect is related to a rapid 
variation of I with the temperature. 

In Ref. magnetic field change experiments have 
been used to determine experimentally the SG coherence 
length in the Heisenberg samples #4 and #5, that was 
compared with numerical simulations corresponding to 
an Ising SG. However, the experimental determination 
of this length scale in the Ising case was not reported 
before. We have measured the zero field cooled magne- 
tization relaxation (mirror procedure of the TRM), after 
an isothermal aging, for various amplitudes of the mag- 
netic field H applied along the c-axis of the Ising sample 
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FIG. 3: Effective waiting times (log scale) derived from the 
field change experiments in the Ising sample #1 as a function 
the magnetic field H . 

1. The effective age t & JJ of the system is evaluated 
by a scaling procedure [12j on reference curves obtained 
with a small field (H = 50e). The effective times as a 
function of the magnetic field were measured for 4 dif- 
ferent waiting times and at two temperatures 0.72T g and 
0.92T g . On a semi- log scale, the data behaves rather lin- 
early with the field (see Fig. 01 . 

The effect of a field variation is to reduce the barri- 
ers that the system overcomes during aging through the 
Zeeman energy E z . The apparent age of the SG, 
smaller than the actual waiting time t w , is then given, 
for small enough E Zl by: 

In t e J f /t w = -E z /k B T. (3) 

The key assumption is that the Zeeman term E z = M • H 
depends on the number N of dynamically correlated spins 
after t w . More precisely, for Ising spins and for not too 
large fields, M should be y/N times the permanent mo- 
ment m/iB of one spin, implying that E z is proportional 
to H and VN (we used a moment of 4(j>b per spins ob- 
tained from the Curie constant above T g ). Following 
Ref. ^3 an d f° r the sake of simplicity, we then simply 
relate the number N of dynamically correlated spins to 
the coherence length through N = a£ 3 . 

The results of the present study along with those repro- 
duced for the Ag:Mn and the thiospinel samples [llj are 
plotted in Fig. 0] (in the Heisenberg case, E z oc N and N 
was extracted in this way in Ref. [llj ) . For a given wait- 
ing time and reduced temperature, the obtained num- 
ber of correlated spins in the Ising sample is noticeably 
smaller than in the Heisenberg samples studied so far, but 
this number grows faster with time, indicating a slower 
growth of the energy barriers with size in the Ising sam- 
ple. The inset of Fig.0]shows the evolution of N with the 
temperature in samples #1 and #4. Clearly, the sepa- 
ration of the active length scales with temperature ( "mi- 
croscope effect" of the temperature @) is much weaker 
in the Ising sample. This is in complete agreement with 
the weak dependence on AT found in the temperature 
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FIG. 4: Number of correlated spins extracted from field 
change experiments in samples #1 at 0.72T 9 and 0.92T 9 , #4 
at 0.72T 9 , 0.78T 9 and 0.84Tg, and #5 at 0.67T 9 , 0.77T 9 and 
0.86T 9 . The curves are simultaneous fits to both field varia- 
tion (Eq. and temperature cycling (Eq. Q and Eq. ex- 
periments, using a single set of parameters for each sample 
(table [IJ. In the main figure, each curve segment is obtained 
at fixed T as a function of t w (in this plot, the universal power 
law of Ref. fTH is a single straight line). The inset shows TV 
as a function of temperature after t w = 1000s for samples #1 
and #4, emphasizing their different behaviors. 
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FIG. 5: The effective waiting times deduced from the tem- 
perature cycle experiments are reproduced from figure^ The 
solid and dashed lines are obtained from Eq. (0 and Eq. Q 
with the very same parameters used in Fig. 0] (see Tabic QJ. 

cycle experiments (small value of dl/d In T in Eq. (J2J). 

Now we compare quantitatively the two sets of exper- 
iments performed in this study. In Ref. |ll|. the cor- 
relation lengths were fitted by a universal power law 
£ = b(t w /T ) aT / Tg (to is the microscopic attempt time 
and is fixed to 10 _12 s), in apparent agreement with nu- 
merical simulations performed on Ising systems [13] . This 
purely activated scenario would correspond to a single 
straight line in Fig.^J and by no means accounts for our 
data in the Ising sample with the parameters used for the 
Heisenberg samples. Moreover, using Eq. one would 
find a slope — ln^/Vo), independent of both a and b, for 
the curve f 2 ff /t 2 vs (AT/T) shown in Fig. □ This sin- 
gle value cannot reproduce the super-activated dynamics 
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TABLE I: Parameters used in Eq. to fit the tempera- 
ture cycle (Fig. Q3 and the magnetic field change experiments 
(Fig.^J. The only free parameters are To, ip, and z since the 
product zv is known (taken from Ref. [To| and Q) and to and 
a have been fixed (see text). The data column gives the total 
number of data points that are fitted for each sample. 
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of the Heisenberg SGs (steep slope in Fig. nor the 
nearly critical dynamics of the Ising sample (weak slope 
in Fig. ^) which requires a renormalized attempt time 
Tq ^> To. Hence, the experimental data have led us to 
propose a crossover expression HE E3 (see also 0): 



^=r^exp (A(T)^/T) 



(4) 



where t' q = tq£ z is a renormalized attempt time (critical 
dynamics) and A(T) = T T g (1 - T/T g f v is the typical 
barrier height that vanishes at T g (super-activated dy- 
namics), following Ref. 7]. As usual, z is the dynamic 
critical exponent and v is the critical exponent that gov- 
erns the equilibrium correlation length. For each of the 
three SG samples #1, #4 and #5, we have fitted the 
whole set of experimental results to Eq. 0] These include 
all field change measurements shown in Fig. 0] and all 
temperature cycling experiments shown in Fig. [3] (in this 
case Eqs.^and 0]are used). A unique set of parameters 
is able to account for all the properties of each sample 
(see table HJ. In order to limit the degrees of freedom in 
our fit, we have imposed the zv values, taking those de- 
rived from dynamic critical scaling (see in refs in Q). We 
have also fixed the geometrical factor a to 2 (N = al 3 ). 
The parameters are not defined with a great quantitative 
accuracy, since their effects on the fit are strongly corre- 
lated. However, a consistent qualitative picture emerges. 
The main tendency is an increasing value of the barrier 
height parameter T and a decreasing barrier exponent ip 
for increasing values of the anisotropy. This behavior of 
the exponent ip is similar to that found in the analysis of 
ac temperature cycling experiments Q, but it contrasts 
with that derived from the scaling of x" relaxations 
as well as with that found in the numerical simulation of 
Ref. ^ij- Note, however, that these previous determina- 
tions were based on a less constrained analysis and, for 
the numerical work, on smaller time scales. 

With a single expression for the growth of the active 
length, Eq. that interpolates between critical and ac- 
tivated dynamics, we have thus been able to reproduce 
consistently two totally different sets of experiments, at 
different temperatures, with a unique set of parameters 
(for each sample). This allows us to estimate the barrier 



height To and the barrier exponent ip, which are no- 
toriously hard to measure. As the anisotropy decreases 
from Ising to Heisenberg SGs, we find an increasingly fast 
separation of the active length scales with temperature, 
corresponding to an increased sharpness of the memory 
effect. Besides, we found that the extracted coherence 
length is noticeably smaller in the Ising sample (large 
To) , but grows faster with time (small ip). At present, 
it is not clear how the strong single spin anisotropy in 
the Ising sample gives rise to both a high value of energy 
barriers and a very small value of the barrier exponent. 
Within a droplet description, ip ~ would imply that 
the droplet energy exponent 9 is also zero, in agreement 
with many recent numerical works on excitations in Ising 
SGs 0|- O ur results underline the role of anisotropy in 
the nature of a SG phase with Heisenberg spins which 
is presently the subject of active investigations (see e.g. 
Ref. HHQ). 
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